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Abstract 



We study the effects of noncommutativity and deformed Heisenberg algebra on the evolution 
of a two dimensional minisuperspace cosmological model in classical and quantum regimes. The 
phase space variables turn out to correspond to the scale factor of a flat FRW model with a positive 
cosmological constant and a dilatonic field with which the action of the model is augmented. The 
exact classical and quantum solutions in commutative and noncommutative cases are presented. 
We also obtain some approximate analytical solutions for the corresponding classical and quantum 
cosmology in the presence of the deformed Heisenberg relations between the phase space variables, 
in the limit where the minisuperspace variables are small. These results are compared with the 
standard commutative and noncommutative cases and similarities and differences of these solutions 
. are discussed. 
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1 Introduction 



As is well known, standard cosmological models based on classical general relativity have no convinc- 
ingly precise answer to the question of the initial conditions from which the universe has evolved. 
This can be traced to the fact that these models suffer from the presence of an initial singularity, 
^ ' the so-called Big Bang singularity. Indeed, there are various forms of singularity theorems in general 

^ . relativity [1], which show that quite reasonable assumptions lead to at least one consequence which 

is physically unacceptable. Any hope of dealing with such singularities would be in the development 
of a concomitant and conducive quantum theory of gravity [2]. On the other hand, one of the most 
important features of theories which deal with quantum gravity is the existence of a minimal length 
below which no other length can be observed [3]. From perturbative string theory point of view, such 
a minimal length, of the order of Planck scale, is due to the fact that strings cannot probe distances 
smaller than the string size. Also, the existence of this minimal length has been suggested in loop 
quantum gravity [4], quantum geometry [5] and black hole physics [6]. Indeed, at the scale of such 

a minimum size i.e. the scales of the order of the Planck length, L = a/^t, the quantum effects of 



gravitation become as important as the electroweak and strong interactions. Clearly, at low energy, 
these quantum gravity effects are not too important, but in high energy physics, that is, energies of 
the order of Planck mass rrip = fi/lp such as the very early universe or in the strong gravitational 
fields of a black hole, one cannot neglect these effects. 

One of the most important features of the existence of a minimal length is that such a length 
is related to what is known as the Generalized Uncertainty Principle (GUP); the usual Heisenberg 
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uncertainty principle should be reformulated at the Planck scale [7, 8]. In a one dimensional system, 
the simplest form of the GUP which shows the appearance of a minimum position uncertainty can 
be written as [7] 



where /? and 7 are positive and independent of Ax and Ap, but may in general depend on the 
expectation values {x) and {p). If we take 7 = it is possible to realize equation (1) from the 

following commutation relation between position and momentum operators 



In [7]- [9], more general GUPs are considered. In more than one dimension GUP naturally implies 
a noncommutative geometric generalization of position space [7]. Noncommutativity between space- 
time coordinates was first introduced by Snyder [10] which has lead to a great deal of interest in this 
area of research in the recent past [11]. 

It is generally an accepted practice to introduce GUP or noncommutativity either through the 
coordinates or fields which may be called geometrical or phase space deformation respectively [12]- 
[17]. Applying GUP or noncommutativity to ordinary quantum field theories where the geometry 
is considered to obey such deformations are interesting since they could provide an eff'ective theory 
bridging the gap between ordinary quantum field theory and string theory, currently considered as 
the most important choice for quantization of gravity. A different approach to GUP and noncommu- 
tativity is through its introduction in the phase space constructed by minisuperspace fields and their 
conjugate momenta [13]- [17]. Since cosmology provides the ground for testing physics at energies 
which are much higher than those on Earth, it seems natural to expect the effects of quantum gravity 
be observed in this context. Alternatively, in cosmological systems, since the scale factor, matter 
fields and their conjugate momenta play the role of dynamical variables of the system, introducing 
GUP and noncommutativity in the corresponding phase space is particularly relevant. 

In general, as we mentioned above, GUP and noncommutativity in their original form (see [7]) 
imply a noncommutative underlying geometry for space-time. However, formulation of gravity in 
a non-commutative space-time is highly nonlinear, rendering the setting up of cosmological models 
difficult. Here, our aim is to study the aspects relating to the application of GUP and noncommuta- 
tivity in the framework of quantum cosmology, i.e. in the context of a minisuperspace reduction of 
dynamics. As is well-known in the minisuperspace approach to quantum cosmology which is based 
on the canonical quantization procedure, one first freezes a large number of degrees of freedom by the 
imposition of symmetries on the spatial metric and then quantizes the remaining ones. Therefore, in 
the absence of a full theory of quantum gravity, quantum cosmology is a quantum mechanical toy 
model (finite degrees of freedom) providing a simple arena for testing ideas and constructions which 
can be introduced in quantum general relativity. In this respect, the GUP approach to quantum 
cosmology appears to be based on physical grounds. In fact, a generalized uncertainty principle can 
be immediately reproduced by deforming the canonical Heisenberg algebra. In other words, the GUP 
scheme relies on a modification of the canonical quantization prescriptions and, in this respect, can be 
reliably applied to any dynamical system. In this sense, one can introduce noncommutativity between 
different dynamical variables of the corresponding minisuperspace and, of course, get different results. 
Here, we rely on and use the most common and accepted practices which have been appearing in the 
literature over the past years. It is to be noted that our presentation does not claim to clear the role 
of GUP and noncommutativity in cosmology in a fundamental way since we study the problem in a 
simple model. However, this may reflect realistic scenarios in similar investigations which deal with 
such problems in a more fundamental way. 

We begin with a flat FRW metric, a positive cosmological constant and a homogeneous scalar 
dilatonic field. We then write the action in the string frame which leads us to a point like Lagrangian 
for the model. We see that the corresponding minisuperspace constructed by the scale factor a and 
dilaton field is curvilinear. Setting up a deformed phase space formalism in such a minisuperspace 
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is not an easy task. Therefore, we introduce a change of variables (a, 0) {u, v) which reduces 
the minisuperspace to a linear (Minkowskian) one. These variables are thus suitable candidates for 
introducing noncommutativity and GUP in the phase space of the problem at hand and enable us 
to present exact solutions for the classical and qTiantum commutative and noncommutative cosmol- 
ogy studied here. Also in the case when the minisuperspace variables obey the GUP commutation 
relations, we obtain approximate analytical solutions for the corresponding classical and quantum 
cosmology. Finally, we compare and contrast these solutions at both classical and quantum levels. 



2 The preliminary setup 

In the pre-big bang scenario, based on the string effective action [18], the birth of the universe is 
described by a transition from the string perturbative vacuum with weak coupling, low curvature and 
cold state to the standard radiation dominated regime, passing through a high curvature and strong 
coupling phase. This transition is made by the kinetic energy term of the dilaton, an scalar field (/) to 
which the Einstein-Hilbert action of general relativity is coupled, see [19] for a more modern review 
and [20] for some exact solutions of string dilaton cosmology. According to this model the lowest 
order gravi-dilaton effective action, in the string frame, can be written as [21] 

1 ^.4. 



^ = -2xJ ^^V^e-^ (7^ + 5^05^0 - 2A) , (3) 

where (l){t) is the dilaton field, is the fundamental string length Ig parameter and A is a (positive) 
cosmological constant. We consider a spatially flat FRW spacetime which is specified by the metric 

= -df + c? {t)6ijdx'dx^ , (4) 

where a{t) is the scale factor. The Ricci scalar corresponding to metric (4) is 

7^ = 6^ + 6^, (5) 

where a dot represents differentiation with respect to t. By substituting (5) into (3) and integrating 
over spatial dimensions, we are led to an effective Lagrangian in the minisuperspace = (a, (p) 

L = e-^ (ed^a - Gda^c^ + a^.^^ - 2Aa^) . (6) 

The momenta conjugate to the dynamical variables are given by 

P„ = _ = e-'^(12da-6aV), P4, = ^ = e''!' {-Qaa^ + 2a^ (f)) , (7) 
oa d(f) 

leading to the following Hamiltonian 

H = \g^^PaPb +U{Q^) = (-^Pf - + ^PaP^) + 2Aa^e-'^- (§) 

Now, it is easy to see that the corresponding minisuperspace has the following minisuper metric 

GABdQ^dQ^ = e-^ i^Uada^ + Ua^dadcj) + 2a^d<i?'^ . (9) 

To apply the deformed commutators to the dynamical variables in a minisuperspace which is repre- 
sented by a curved manifold with a minisuper metric given by (9), in a natural generalization, one can 
replace in (2) with Q^^PaPb- In general, this generalization does not provide a suitable expression 
because of the ambiguity in the ordering of factors Q and Pq. Therefore, the above minisuperspace 
does not have the desired form for introducing noncommutativity and GUP among its coordinates. 



3 



To avoid the physical difficulties and simplify the model, consider the following change of variables 

= (a, (P)^q^ = (u, v) 

u + v = Aa'^e~^'^, u - v = a^e''^''^ , (10) 
where a and /3 are two constants which satisfy the relations 

3 

a + /3 = 3, a/? = -. 
In terms of these new variables, Lagrangian (6) takes the form 

with the corresponding Hamiltonian becoming 

n = \{pl-pt)+u^{u^-v^), (12) 

which describes an isotropic oscillator-ghost-oscillator system with frequency uP' = ^. Thus, in the 
minisuperspace constructed by = {u, v) , the metric is Minkowskian and represented by 

GABdq^dq^ = ^ (du^ - dv^) . (13) 

Now, we have a set of variables (u, v) endowing the minisuperspace with a Minkowskian metric and 

hence this set of dynamical variables are suitable candidates for introducing noncommutativity and 
GUP in the phase space of the problem at hand. The final remark about the above analysis is that 
Lagrangian (6) possesses an interesting symmetry thanks to the presence of the stringy dilaton. This 
symmetry exhibits itself through the transformation [21] 

a{t) l/a(i), (t){t) (pit) - 61na(i). (14) 

It is easy to show that Lagrangian (6) is invariant under this transformation. Such symmetry (duality) 
is one of the major features of the solutions of equations of motion in string dilaton cosmology [20], 
so that if the set of variables (o, 0) solve the equations of motion, the set (1/a, (/> — 6 In a) is also a 
solution. On the other hand, in terms of the variables (-u, v) the Lagrangian takes the simple form 
(11) yielding linear differential equations for the corresponding dynamical equations. Therefore, the 
duality symmetry is nothing but a suitable linear combination of u and v. Indeed, one can easily show 
that Lagrangian (11) and also Hamiltonian (12) are invariant under the following transformations ^ 

17 15 15 17 

—u - —V, —u - —V, (15) 

17 15 15 17 

Pu^ -^Pu + yP-"^ Pv ^ --^Pu- -^Pv (16) 

The preliminary setup for describing the model is now complete. In what follows we will study 
the classical and quantum cosmology of the minisuperspace model described by Hamiltonian (12) in 
noncommutative and GUP frameworks. 



^In general, Lagrangian (11) is invariant under pseudo rotations in two dimensional Minkowskian space 

u — » w cosh § -\- V sinh §, v —>■ u sinh + v cosh 

where is the parameter of transformations. In (15) and (16) we take a special choice for '& to recover the duality of 
the theory represented by (14). 
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3 Classical model 



As mentioned above, the dynamical system described by Hamiltonian (12) is a simple isotropic 
oscillator-ghost-oscillator system and its classical and quantum solutions can be easily obtained. 
Since our aim here is to study the effects of deformed Poisson brackets on the classical trajectories, in 
what follows we consider commutative, noncommutative, and GUP classical cosmologies and compare 
the results with each other. In the next section we shall deal with the quantum cosmology of the 
model. 

3.1 Classical cosmology with ordinary Poisson brackets 

As is well known for a dynamical system with phase space variables {qi,Pi), the Poisson algebra is 
described by the following Poisson brackets 

{Qi, Qj} = {Pi,Pj} = 0, {qi,Pj} = Sij, (17) 

where in our case qi{i = 1,2) Therefore, the equations of motion 

become ^ 

u = {u,n} = -pu, pu = {Pu^T-L} = -2u?u, (18) 

v = {v,H} = ~Pv, Pv = {py,H} = 2u^v. (19) 



Integrating the above equations, one is led to 



u{t) =lul + ^\ sin 



1/2 

-1 ^Uq 

Lot + tan 



Uq 



v{t) = + ^ sin 



1/2 

ut + tan 



where for the initial conditions we take 



u{t = 0) = uo, u{t = 0) = iiQ, 
v{t = 0) = vo, v{t = 0) = Vq. 



(20) 
(21) 

(22) 



The above solutions must satisfy the Hamiltonian constraint, H = 0. Thus, substitution of equations 
(20) and (21) into (12) gives the following relation between integration constants 

U0 + ^=V0 + ^- (23) 
Prom the above equations, we see that the classical trajectories obey the relation 

. = ±cosftan-i^-tan-i^^V±s-ftan-i — -tan-^— ) l^ + ^.-^'Y^- (24) 
V ^^0 uo J \ Vo no / \ y 

Note that the minisuperspace of the above model is a two-dimensional manifold which in terms of the 
old variables a and (p is represented by < a < oo, — oo < (p < +oo. Following [22], we may divide its 
boundary into two, the nonsingular and singular. The nonsingular boundary is the line a = with 
1^1 < +00, while at the singular boundary, at least one of the two variables is infinite. In terms of the 
variables u and v, introduced in (10), the minisuperspace is recovered by li > 0, —u < v < u, and the 
nonsingular boundary may be represented hy u = v = 0. This discussion leads us to the imposition 
of more restrictions on the initial conditions (22) such that the classical trajectories would no longer 
meet the nonsingular boundary. This condition is achieved when the coefficient of the second term 
in (24) is nonzero. 
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Now, let us go back to the old variables a and (p, in terms of which we obtain the corresponding 
classical cosmology as 

4a(tr-^ = 4a(t)±^=^, 

u — V 

(j){t) = 2(3ln\a{t)\ -2ln\u-v\, 
leading to the following sets of classical solutions 



(25) 
(26) 



and 



a+(t) = ao [tana;(i — to)]^^^ , 
4>-\-{t) = {^/3 — 1) In I tanu;(i — io)| — 2 In | cosa;(i — to)\ + (po, 

a-{t) = [aotana;(t - to)]"^^"^, 



(27) 



(28) 



^_(t) = {-VS- l)ln|tana;(t-io)| - 21n | cosa;(i - io)| + (po, 

where ao, to and (f)o are some constants which can be written in terms of uq, uq, vq and vq. These 
two sets of solutions are related to the duality symmetry (14), indeed we have 

1 



a_(t) 



0_(t) = (j)+{t) - 61na+(t). 



(29) 



As we mentioned in the previous section, in the minisuperspace (u, v) the duality symmetry is denoted 
by the linear combination (15) of u and v. Therefore, applying the duality transformation (15) we 
are led to the following set of solutions 



1/2 



U{t)=lut, + ^ -sin 



V,2 



17 



1/2 



Vit) 



^0 + ^ I S — sm 



15 



Lot + tan 



ojt + tan 





15 


uo . 




-1 i^uo' 


17 







■ sm 



Lot + tan 



Lot + tan" 



-1 ^ 



(30) 

(31) 



It is clear that these solutions are essentially a special linear combination of (20) and (21) which 
obviously solve the classical equations of motion because of their linearity. 



3.2 Classical cosmology with noncommutative phase space variables 

Let us now proceed to study the behavior of the above model in a deformed phase space framework 
such that the minisuperspace variables do not (Poisson) commute with each other. In general, non- 

commutativity between phase space variables can be understood by replacing the usual product with 
the star-product, also known as the Moyal product law between two arbitrary functions of position 
and momentum as 

f{xi)g{x2)\ (32) 

where a"'^ denote the noncommutative parameters [23]. Here, we consider a noncommutative phase 
space in which the Poisson algebra is a deformed one given by 

{QinciQjnc\ ~ ^^iji {PinctPjnc} ~ 0, {Qinci Pjnc} ~ ^iji (^^) 

where e^j and 5ij are Levi-Civita and Kronecker symbols respectively, qinc{i = 1,2) = Unc^Vnc and 
Pinci} = 1)2) = PunciPvnc- With the deformed phase space defined above, one may consider the 
Hamiltonian of the noncommutative model as having the same functional form as (12), but with the 
dynamical variables satisfying the above-deformed Poisson brackets, that is 

^nc = \ {pI„, - pL) + («L - vL) ■ (34) 



(/ *a 9){x) =exp 
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Thus, the dynamics of the system can be described by the following equations of motion 



U'nc = {Unc, T-(-nc} = ^Pur,c " "^Ooj^nc, Punc = {Punc^Unc} = -2uPunc, (35) 
V'nc = {Vnc, Ti-nc} = - ]^Pvuc " ^Ouj'^Unc, Pvnc = {Pvnc^'^nc} = '^'^'^Vnc- (36) 

Eliminating the momenta from the above equations, wc get 

Unc + ^"^Unc + ^Ou^Vnc = 0, (37) 
V'nc + i^'^Vnc + ^OuP'Unc = 0. (38) 

We see that the noncommutative parameter appears as a coupling constant between equations of 
motion for Unc and Vnc- Integrating equations (37) and (38) yields 

Unc{t) = ^e^"''* sin [to^l - O^uH + 5i\ + ^e"^'^"* sin [uj^/l - B^ujH + 82] , (39) 

Vnc{t) = -^e^'^'* sin [wVl - O^ioH + Si] + He"^'^"* sin [uVl - e^uH + 62] , (40) 

where A, B, Si and 62 are integrating constants. The requirement that the noncommutative Hamil- 
tonian constraints should hold during the evolution of the system, that is, Ti-nc = 0; leads to the 
following relation between integrating constants 

^ = 0. (41) 

This means that in the noncommutative minisuperspace the system follows one of the trajectories 
Unc = Vnc (if .4 = 0) Or Unc = —^nc (if B=0). In the case when ^ = 0, the two coordinates behave 
similar to two coupled springs, oscillating back and forth together like and f— with frequency 
wVl — and an exponentially damping amplitude. Alternatively, if we take B = 0, the two 

variables oscillate in opposite directions like < > and — ><— with the same frequency wVl — O'^oo'^ but 

with an exponentially increasing amplitude. 

Before going any further some remarks are in order. An important ingredient in any model theory 
related to the quantization of a cosmological setting is the choice of the quantization procedure used to 
quantized the system. The most widely used method has traditionally been the canonical quantization 
method based on the Wheeler-DeWitt (WD) equation which is nothing but the application of the 
Hamiltonian constraint to the wavefunction of the universe. A particularly interesting but rarely used 
approach to study the quantum effects is to introduce a deformation in the phase space of the system. 
It is believed that such a deformation of phase space is an equivalent path to quantization, in par 
with other methods, namely canonical and path integral quantization [24]. This method is based on 
Wigncr quasi-distribution function and Weyl correspondence between quantum mechanical operators 
in Hilbert space and ordinary c-number functions in phase space. The deformation in the usual phase 
space structure is introduced by Moyal brackets which are based on the Moyal product (32) [23]. 
However, to introduce such deformations it is more convenient to work with Poisson brackets rather 
than Moyal brackets. 

Prom a cosmological point of view, models are built in a minisuperspace. It is therefore safe to 
say that studying such a space in the presence of deformations mentioned above can be interpreted as 
studying the quantum effects on cosmological solutions. One should note that in gravity the effects of 
quantization are woven into the existence of a fundamental length [3] . The question then arises as to 
what form of deformations in phase space is appropriate for studying quantum effects in a cosmological 
model? Studies in noncommutative geometry [11] and GUP [7] have been a source of inspiration 
for those who have been seeking an answer to the above question. More precisely, introduction of 
modifications in the structure of geometry in the way of noncommutativity has become the basis from 
which similar modifications in phase space have been inspired. In this approach, the fields and their 
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conjugate momenta play the role of coordinate basis in noncommutative geometry [25]. In doing so 
an effective model is constructed whose validity will depend on its power of prediction. 

A question worth asking at this stage is: would the noncommutative scheme presented above 
really offer a quantum picture of the model at hand and should we refrain from using any other 
quantization method simultaneously? It is important to note that equivalence between the two 
different approaches of quantization cannot hold true in models where deformation in phase space 
is introduced in a Lorentz non-invariant manner, like what we have done here. This is not hard 
to understand since the WD equation is a direct consequence of diffcomorphism invariance and so 
if a deformation in phase space breaks such an invariance then the results of different quantization 
methods should be different. For models where the Lorentz invariance deformation is studied see [26]. 
Therefore, in the next section when we quantize our model we also invoke noncommutative quantum 
cosmology based on the star-product WD equation. 

3.3 Classical cosmology vi^ith GUP 

In more than one dimension a natural generalization of equation (2) is defined by the following 
commutation relations [7] 

[xi,pj] = i (Sij + pSijp^ + iS'piPj'j , (42) 
where = J^PiPi P,P' > are considered as small quantities of first order. Also, assuming that 

\pi,Pj]=0, (43) 
the commutation relations for the coordinates are obtained as 

.(2/3 -/?') + (2/? + /3')/3p' 



{piXj - pjXi) . (44) 



As it is clear from the above expression, the coordinates do not commute. This means that to 
construct the Hilbert space representations one cannot work in the position space. It is therefore more 
convenient to work in momentum space. However, since in quantum cosmology the wavefunction of 
the universe in momentum space has no suitable interpretation, we restrict ourselves to the special 
case P' = 2/3. As one can see immediately from equation (44), the coordinates commute to first order 
in /? and thus a coordinate representation can be defined. Now, it is easy to show that the following 
representation of the momentum operator in position space satisfies relations (42) and (43) (with 
/?' = 2/?) to first order in /3 

Pi = -z[l- -— ^ — . (45) 



Equations (42)- (45) may now be realized from the following commutation relations between position 
and momentum operators 

[u, pu] = i{l + pp^ + 2PpI) , [v,py]=i{l + pp^ + 2PpI) , (46) 

[u,Pv] = [v,Pu\ = ^i(3PuPv, (47) 

[xi,Xj] = \pi,Pj] = 0, Xi{i = 1,2) = u,v, pi{i = 1,2) =Pu,Pv (48) 

Before quantizing the model within the GUP framework in the next section, we investigate the effects 
of the classical version of GUP, i.e. the classical version of commutation relations (46)- (48) on the 
above cosmology. As is well known, in the classical limit the quantum mechanical commutators should 
be replaced by the classical Poisson brackets as [P,Q] — > ih{P,Q}. Thus, in classical phase space 
the GUP changes the Poisson algebra (17) according to 

{u,p^} = 1 + /?/ + 2l3pl {v,p^} = 1 + /3p2 + 2/3pl, (49) 
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Figure 1: Approximate behavior of u{t) and v{t) in the classical GUP framework. The figures are plotted for rmmerical 
values uj = 1 and /3 = 0.01. We take the initial conditions u{t = 0) = v{t = 0) = 0, pu{t = 0) = pv{t = 0) = 3 for the 
figure on the left and u{t = 0) = v{t = 0) = 0, pu{t = 0) = —pv{t = 0) = 3 for the figure on the right. 



{u,Pv} = {v,Pu} = 2/3puPv, (50) 

{xi,Xj} = {pi,pj} = 0, Xi{i = 1,2) = u,v, pi{i = 1,2) =pu,Pv, (51) 

where p^ = ^{Pu~ pD- Such deformed Poisson algebra is used in [27] to investigate the effects of 
deformations on the classical orbits of particles in a central force field and on the Kepler third law. 
Also, the stability of planetary circular orbits in the framework of such deformed Poisson brackets is 
considered in [28]. Note that here we deal with modifications of a classical cosmology that become 
important only at the Planck scale where the classical description is no longer appropriate and a 
quantum model is required. However, before quantizing the model we shall provide a deformed 
classical cosmology. In this classical description of the universe in transition from commutation 
relation (2) to its Poisson bracket counterpart we keep the parameter /3 fixed as ^ — 0. In string 
theory this means that the string momentum scale is fixed when its length scale approaches zero. 
Therefore, the equations of motion read 

u = {u, n} = ^pu (l + 5l3p^) , Pu = {Pu, n} = -2uj^u (l + fip" + 2ppl) + ^uj'^PvpuPv, (52) 

v = {v,n} = - ^py (l - 3/?p2) , = n} = 2j'v (l + (ip'' + 2ppl) - Au^PupuPv (53) 

We see that the deformed classical cosmology forms a system of nonlinear coupled differential equa- 
tions which unfortunately cannot be solved analytically. In figure 1, employing numerical methods, 
we have shown the approximate behavior of u{t) and v{t) for typical values of the parameters and 
initial conditions. As it is clear from the figure, both variables u and v repeat their back and forth 
oscillatory behavior as the noncommutative case, but here the oscillations are not harmonic. Like 
the noncommutative case, depending on the initial conditions, the minisuperspace variables behave 

as and <— <— , (see the right figure) or < > and ^•f— (see the left figure) i.e. they move back and 

forth either in the same or in opposite directions. A comment on the above results is that although 
in the limit /3 ^ one can recover the ordinary classical cosmology described by equations (18) 
and (19), as this figure shows, taking a nonzero value for /? may disturb the oscillatory nature of 
the universe. Also, in the presence of /? terms, the period of oscillations becomes larger and thus 
the Big-Crunch in the corresponding cosmological model occurs later in comparison with the usual 
models where (3 = 0. This means that the effects of GUP are important not only in the early but also 
at late times in the cosmic evolution. In fact, within the GUP framework, the quantum gravitational 
effects may be detected at large scales as well. 

4 The quantum model 

We now focus attention on the study of the quantum cosmology of the model described above. Here, 
as in classical cosmology, for comparison purposes between ordinary commutative, noncommutative. 
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and GUP frameworks, we consider the quantum cosmology of the model separately in each case and 
compare the results. Our starting point is to construct the WD equation from the corresponding 
Hamiltonian. 



4.1 Quantum cosmology with ordinary commutation relations 

The quantum version of the model described by relations (17) can be achieved via the canonical 
quantization procedure which leads to the WD equation, 7^^' = 0. Here, H. is the operator form of 
the Hamiltonian (12) which annihilates the wavefunction ^. By replacing pq — > —i-§^ in (12) we get 
the WD equation as 



^{u,v)=Q. (54) 



This equation is a quantum isotropic oscillator-ghost-oscillator system with zero energy. Therefore, 
its solutions belong to a subspace of the Hilbert space spanned by separable cigenfunctions of a two 
dimensional isotropic simple harmonic oscillator Hamiltonian. Separating the cigenfunctions of (54) 
in the form 

yields 

UnM = [-) -^^—^H^A^ul (56) 

K.W=(-) ^=^H^.iV2^vl (57) 

subject to the restriction ni = n2 = n. In (56) and (57), Hn{x) are Hermite polynomials and the 
cigenfunctions are normalized according to 



/. 



+ 0O „ 



Hn{x)Hm{x)dx = T'K^I'^nXb^n- (58) 



Now, we impose the boundary condition on these solutions such that at the nonsingular boundary 
(at n = = 0) the wavefunction vanishes [29], that is, *(n = 0, u = 0) = 0, which yields 

Rn{u = 0)Hn{v = 0) = 0^n = odd. (59) 

In general, one of the most important features in quantum cosmology is the recovery of classical 
cosmology from the corresponding quantum model, or in other words, how can the WD wavefunctions 
predict a classical universe. In this approach, one usually constructs a coherent wavepacket with 
good asymptotic behavior in the minisuperspace, peaking in the vicinity of the classical trajectory. 
Therefore, we may now write the general solution of the WD equation as a superposition of the above 
cigenfunctions 

1 /2 

^{u,v)=(—) e-("'+''') ^,Hn{^u)H^{^v). (60) 

n=oda 

Figure 2 shows the square of the wavefunction and its contour plot. As we can see from this figure 
the peaks follow a path which can be interpreted as the classical trajectories (24). The crests are 
symmetrically distributed around v = which correspond to the it signs in (24). Thus, it is seen 
that there is an almost good correlations between the quantum patterns and classical trajectories in 
the u — V plane. 
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0.5 1 1.5 2 



Fi gure 2: The figure on the left shows |\['(u, the square of the commutative wavefunction while the figure on the 
right, the contour plot of \'^{u, The figures are plotted for numerical value a; = 1 and we have taken a superposition 
of eight terms in (60) with all c„ up to cis taken to be unity. 



4.2 Noncommutative quantum model 

Now, the quantum version of the noncommutative cosmology is achieved by replacing the Poisson 
brackets with the corresponding Dirac commutators, {,} — > —*[,]• Thus, the Poisson brackets (33) 
between minisuperspace variables should be modified as follows 

[UncVnc] = iO, [UncPu] = [VncPv] = «■ (61) 

The corresponding quantum cosmology can be obtained by modification of the operator product 
in the WD equation TC^ = with the Moyal deformed product Tine * ^' = 0, where Tine is the 
noncommutative Hamiltonian (34). However, there is an alternative expression for the Moyal star 
product which is given by the shift formula 

Unc = U- -6py, Vnc = V + -Opu, (62) 
Pn„, =P«, Pv^^=Pv (63) 

It can easily be checked that if the noncommutative variables obey relations (61), then {u,v,pu,Pv) 
satisfy the usual Heisenberg algebra 

[u,v] = 0, [u,pu] = [v,Pv] = i, [u,Pv] = [v,Pu] = 0. (64) 

In terms of these commutative variables the Moyal WD equation Tine * ^ = transforms to the usual 
WD equation TC^ = with Hamiltonian 

n = ^{l- Lo^e^^ (pI - pI) + Lo^ (u^ - - Olo^ {up, + vpu) . (65) 

Therefore, the noncommutative version of the WD equation can be written as 
1 

where 



n^{u,v) = H^{u,v) 



\ ifu - pI) + {u' - v') - en^ {up, + vpu] 



^{u,v) = 0, (66) 
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The Hamiltonian operator consists of two parts 

H = — 9Q^Luv, (68) 

where 



uv 



^ {pI - pI) + (u^ - v"^) , and L^y = up^ + vpu- (69) 
It is easy to see that these two operators commute with each other 

[H^v,Luv] = 0, (70) 

which means that H^v and L^v have simultaneous eigcnfunctions ^. In the coordinate system {u,v) 
the WD equation (66) is not separable. Thus, to transform it to a separable equation, consider the 
new variables (p, (p) defined as 

u = pcosh(^, V = psinhip, (71) 

in terms of which we have 

■9 . d . d 

Luv = -m- IV— = -t^, (72) 

OV OU Oip 

^- = -4(9^-9^]+^ (^ -^ ) = -4 + ^''^ 

Therefore, if we separate the solution of the WD equation according to 5'(p, ip) = R{p)e'^'^^, it is easy 
to see that this is an eigenfunction of Luv with eigenvalue u. The requirement that ^{p,ip) should 
solve the WD equation yields the following differential equation for R{p) 

^ + + ( 4.0^0 - An'p' + ^]r = 0. (74) 
ap^ p dp y p'^ J 

This equation, after a change of variable r = 2Q.p^ and transformation R = ^W, becomes 

(fW ( 1 K 1/4 

where k = ^i'Q,9 and /x = z|. The above equation is the Whittaker differential equation and its 
solutions can be written in terms of confluent hypergeometric functions M(a, 6; x) and U{a, b; x) as 



W{r) = 6-'/^"+^ 



cU{ii- K + ^,2ii + l;r) + c'M{ii - k + ^,2p + l;r) 



(76) 



In view of the asymptotic behavior of M{a,b;x) ~ e^/x^ " [30], we take c' = 0. Thus, going back, 
the eigcnfunctions of the WD equation can be written in terms of variables (p, ip) as 

Mp, ^) = P'^e-'^'^U (1:1^ + z^, 1 + zz.; 2Jlp2) e^^. (77) 



^Like the theory of ordinary harmonic oscillator, we can also define the creation and annihilation operators 



and similarly for dv and dj, satisfying the following commutation relations 

[a„,al] = [at,,al] = 1, 

with other commutators being zero. In terms of these operators H-av and L^v can be viewed as 
Use of the same commutators as above between d and d^ makes it easy to see that [Huv,Luv] = 0. 
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Figure 3: Up: the figure on the left shows v)\'^ , the square of the commutative wavefunction in (p, (p) coordinates 

while the figure on the right shows its contour plot. Bottom: the same figures in the noncommutative case. The 
figures are plotted for the numerical value fl = 1 and we have taken a superposition of three terms in (78) with 
Ci = C2 = Cs = 1. 

If we demand that ^'j/ should be a single-valued function of ip, then 1/ should be an integer v = n. 
Now, we may write the noncommutative wavefunction in the (p, (p) coordinates as 

*(p,^)= £ C7„p-e-^pV(i::|^ + i^,l + m;2^^p2^e^n^. (78) 

n=— 00 ^ ^ 

Note that in terms of variables {p^'-p), the boundary condition ^{u = 0,i; = 0) =0 takes the form 
^(p = 0, tp) =0 which for (78) is automatically held. 

In figure 3 we have plotted the wavefunction (78) in the noncommutative quantum model. In the 
case where = (commutative model), although we have analyzed the behavior of the wavefunction in 
the previous subsection, we have also considered it again in the (p, ip) coordinates. As we can see from 
this figure, the commutative wavefunction has two dominant peaks in the vicinity of p = which then 
follow a path that can be interpreted as the classical trajectory (compare with the similar behavior 
in figure 2). Therefore, like the commutative wavefunction in the {u,v) coordinates, there is also a 
correlation between the classical and quantum schemes. On the other hand, the noncommutative 
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wavefunction predicts the emergence of the universe from a state corresponding to one of the two 
dominant peaks. Although there are some small peaks in this figure, as p grows, their amplitude are 
suppressed. We see that the correlation with classical trajectories is missed, i.e. the noncommutativity 
implies that the universe escapes the classical trajectories and approaches a stationary state. 



4.3 Quantum model with GUP 

The study of quantum cosmology of the model presented above in the GUP framework is the goal we 
shall pursue in this subsection. The Hamiltonian of the model is given by (12) and the corresponding 
commutation relations between dynamical variables by equations (46)- (48). In the WD equation we 
take the representation (45) for the momenta Pu and py which leads to the following equation up to 
first order in (3 



We again separate the solutions in the form ^{u,v) = U{u)V{v), leading to 



*(tx,t;)=0. (79) 



where Wi{i = 1,2) = U,V, Wi{i = 1,2) = u,v and A is the separation constant. The appearance of 
a fourth order differential equation to describe a physical phenomena is interesting, since it requires 
investigation of the corresponding boundary conditions which is not the goal of our study in this 
paper. These equations cannot be solved analytically, but since /3 is a small parameter which appears 
only in the fourth order term, we may look for an approximation method which, within its domain 
of validity, leads us to a second order differential equation. From figure 2, we note that the dominant 
peaks of the wavefunction occur in the vicinity u,v ~ 0. On the other hand, the effects of (3 
are important at the Planck scales, or in cosmology language in the very early times of the cosmic 
evolution, which in our model means u,v ~ 0. When /? = 0, the solutions of equation (80) for 
A/2u; = 2ri + 1 are given by (56) and (57). In the limit — > we may take e~'^^" +^ ) ^ l which 
means that Wi oc Hn{V2ojWi). Therefore, in this limit, we have 



d'^Wi /2w\i/^64w2n(n-l)(n-2)(n-3)^^ . ^ 
dwf V vr y V2"n! 

where we have used the relation dHn{x)/dx = 2nHn-i{x). Now, using the following series formula 
for Hermite polynomials 

Hn{x) = E(-l)^ (,_2,s)!.! ('"^""'^' 

we are led to the approximate formula Hn-A{x) oc Hn{x) in the case of a small argument and n > 4. 
Thus, in our approximation, equation (80) becomes 



d^Wi 



+ 2a; (2n + 1 - 2/3o - 2a;K;f ) Wi = 0, (83) 



dwf 

where Pq, up to a numerical factor, is 

/3o = y/?a;n(n - l)(n - 2)(n - 3). 
The solutions of the above equation up to a normalization factor can be written as 

Wn{w) = e-'^"''if„_;3o(v^y;). (84) 



14 



Figure 4: The figure on the left shows v)f, the square of the GUP wavefunction while the figure on the right, the 

contour plot of The figures are plotted for numerical value a; = 1, /9o = 4 and we have taken a superposition 

of four terms in (85) with all Cn taken to be unity. 

Finally, the general solutions of the WD equation (79) for small values of u and v are as follows 



where to recover solution (60) in the limit Pq — > 0, the summation is taken over the odd values of n. 
Also, to satisfy the boundary condition ^'(n = 0, w = 0) = 0, we choose /?o such that n — /3o is an odd 
integer for n > 4. 

Figure 4 shows the wavefunction of the corresponding universe when the minisuperspace variables 
obey the GUP relations for small values of u and v. As is clear from this figure the wavefunction has 
two single peaks which are symmetrically distributed around v = 0. Compare to the commutative 
wavefunction (see figure 2), here we have no wave packet with peaks following the classical trajectories. 
We see that instead of a series of peaks in the ordinary WD approach, we have only a couple of 
dominant peaks. This means that, similar to the noncommutative case and within the context of the 
GUP framework, the wavefunction also shows a stationary behavior. One may then conclude that 
from the point of view adopted here, noncommutativity and GUP may have close relations with each 
other. 

We should note that the above analysis on the behavior of the GUP wavefunction is achieved 
in the region u,v ~ 0. If we relax this approximation, equation (81) is no longer valid and in 
the perturbation method we should use solutions (56) and (57) in the /3-term of equation (80) in 
order to obtain the GUP wavefunction. However, since for large values of u and v these solutions 
have an exponentially decreasing behavior, we do not expect a major effect on the behavior of the 
above wavefunction even if we consider the problem more rigorously. We can also rely on the above 
approximate GUP solutions for a wider range of u and v. 

To clarify this point we investigate the quantum GUP model in a different representation. To do 
this, we remember that the WD equation (79) is based on the representation (45) of momenta in the 
u — V space. In this sense we note that the existence of a minimal length means that we cannot have 
localized physical states. However, as is shown in [7], when there is no minimal uncertainty in mo- 
mentum one can work with the momentum space wavefunction ^{p) with the following representation 



(85) 



n=odd 



[7] 



(86) 
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(87) 



where in our model we have, as before Xi{i = 1,2)= u, v, pi{i = 1,2)= Pu,Pv and = \{p\ — pD- 
Thus, one can now define the proper physical states with maximal localization and use them to define 
a "quasi-position wavefunction" . In [7] it is shown that the transition between the momentum space 
representation of the wavefunction and its quasi-position counterpart is given by a generalized Fourier 
transformation as 



dp 



ixtan" 



(1 + /?p2)3/2 



exp 



VP J 



Hp)- 



(88) 



Now, using representations (86) and (87), we get the following form for the WD equation in momentum 
space 



d d \ p^ 

^"'dfu^^'^Wv) ^ 2a;2(l + ^p2)2 



^ijPu,Pv)=^, 



(89) 



where in contrast to the minisuperspace coordinate representation of the WD equation, it is a second 
order differential equation. This equation, up to first order in (i and also neglecting in the last 
term takes the form ^ 



5^ 

dpi dpi 



d'^ d d 

V OPu OPy 



1 



+ ^{pI-pI) 



HPu,Pv) = 0, (90) 
which is a separable equation and its solutions may be written as ^{pu,Pv) = n(p„)T(p„), leading to 



(f ^ d f 1 
+ PPu ' — 



dpt 



dpi 



- Ppv 



dpv 
_d_ 

dpy \ Alu^ 



f ' 



^{Pu) = 0, 



Pv + E 



T{pv) 



0, 



(91) 



(92) 



where £^ is a separation constant. The above equations have exact solutions in terms of confluent 
hypergeometric functions as 



^e{Pu) =Pwexp 



4w 



P] 2 
T 



4oj 



U 



U 



Pv , 



4 A^TTW^ 2 



2u 



(93) 
(94) 



Therefore, up to first order in /3 the eigenfunctions of equation (90) can be written as 



^E{Pu,Pv)=PuPve U. + 4)f«e 4)P^ul-+ ^ '2'£) li 



3 uA2E-Jl 3 pi_ 
^ 4 '2'2o; 



Now the eigenfunctions in the u — v representation can be obtained from (88) 



(95) 



,T/ ( \ r"^ ' 

TT J-oo J-oo (1 



dpudpv 



(l+/3p2)3 



exp 



ill tan ^{^fPPu) 

. TP . 



exp 



ivtan ^{VPPv) 

TP 



^EijPu,Pv)- 

(96) 



^Our classical analysis (see equations (20)-(23), (39)-(41) and figure 1), shows that the dynamical variables u and 

V have the same order of magnitude. Therefore the Hamiltonian constraint "H = Q implies that p„ and p„ should also 
have the same order of magnitude. In this sense, neglecting in equation (89) is quite reasonable. 
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Figure 5: The figure shows the square of the GUP wavefunction (97). The figure is plotted for numerical 

values u) = 1/24, /3 = 8 and we have taken a superposition of four terms in (97) with all C(E) taken to be unity. 



It is seen that this expression is too comphcated for extracting an analytical closed form for the 
eigenfunctions in terms of the minisuperspace variables u and v. Moreover, the wavefunction is a 
superposition of these eigenfunctions as 



On the other hand to construct such a superposition we need to know whether the separating constant 
E has a quantized or a continuum spectrum. Here, we do not intend to deal with such questions 
since they are out of the scope of our study. However, in figure 5 we have plotted the square of the 
wavefunction where we have taken a discrete superposition of four terms in (97), all to first order in 
(3. It should be noted that we have relaxed the assumption of small values for u and v in this analysis 
and figure 5 shows the square of the wavefunction for a wide range of these variables. Although 
this figure does not completely coincide with figure 4, as is clear, the wavefunction has its dominant 
amplitude in regions represented by the small values of u and v, in agreement with our previous 
analysis of the GUP quantum model based on equation (81). 

5 Conclusions 

In this paper we have studied the effects of deformations (noncommutativity and GUP) in phase space 
on the cosmic evolution of a two dimensional minisuperspace model. Our starting point was the lowest 
order gravi-dilaton effective action in the string frame where the Einstein-Hilbert action with a positive 
cosmological constant is augmented with an scalar field, the dilaton. By considering a flat FRW metric 
for the space time geometry in this action, we obtained the corresponding effective Hamiltonian in the 
minisuperspace constructed by the scale factor a and dilaton fleld We saw that this minisuperspace 
has a curved metric and thus does not have the desired form for introducing noncommutativity and 
GUP between its coordinates. For this reason, we introduced a new set of variables (u,^), in terms 
of which the minisuper metric took a Minkowskian form and the Hamiltonian of the model described 
a simple isotropic oscillator-ghost-oscillator system in which the cosmological constant plays the role 
of frequency. Although this is a simple toy model, these variables are suitable candidates for a 
phenomenological study of noncommutativity and GUP in the corresponding phase space. Another 
feature of the model in the (ti, v) coordinates is that the duality symmetry of the string dilaton 
action exhibits itself as a special linear combination of u and v. In the case of (Poisson) commutative 
phase space we saw that both dynamical variables have oscillatory behavior with the same amplitude. 




(97) 



E 
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Depending on the initial conditions, these osciUations may occur in the same or opposite directions. 
As for the quantum version of this commutative model, we obtained exact solutions of the WD 
equation. The wavefunction of the corresponding universe consists of two branches where each may 
be interpreted as part of the classical trajectory. We saw that since the peaks of the wavefunction 
follow the classical trajectory, there seems to be good correlations between the corresponding classical 
and quantum cosmology. 

We also studied the noncommutative cosmology where the commutator between the minisuper- 
space variables is deformed through a noncommutative parameter 9. In the classical noncommutative 
model this parameter plays the role of a coupling constant between equations of motion for dynamical 
variables u and v. We solved the equations of motion exactly and showed that the variables u and 
V oscillate in the same or opposite directions with an exponentially damping or increasing amplitude 
respectively. In the noncommutative quantum cosmology, we found that the corresponding Hamil- 
tonian is constructed out of two operators which commute with each other and thus have the same 
eigenf unctions. Although the ensuing WD equation in this case is not separable in the {u, v) coordi- 
nates we introduced a new set of variables {p, (p) in terms of which the WD equation was amenable 
to exact solutions in terms of confluent hypergeometric functions. Finally, when the phase space 
variables obey the GUP relations wc constructed the classical equations of motion and seen that they 
form a system of nonlinear differential equations which cannot be solved analytically. We showed 
that the oscillatory behavior of u and v in the same or opposite directions is again repeated but in 
this case the oscillations are not harmonic. We also found that the GUP parameter (3 causes a larger 
period in comparison to the ordinary model where /3 = 0. This can be interpreted as the importance 
of quantum gravitational effects not only at early times but also at late times of the cosmic evolution. 
The resulting quantum cosmology and the corresponding WD equation in the GUP framework were 
also studied and approximate analytical expressions for the wavefunctions of the universe were pre- 
sented in the limit of small u and v variables. These solutions show only one possible state with no 
classical correlation at early times from which our universe can emerge. We saw that such behavior 
also occurs in the noncommutative quantum model, showing that from the point of view adopted 
here, noncommutativity and GUP may be considered as similar concepts. 
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